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ABSTRACT 

Pulsar dynamic spectra sometimes show organised interference patterns; these patterns have 
been shown to have power spectra which often take the form of paraboHc arcs, or sequences of 
inverted paraboHc arclets whose apexes themselves follow a parabolic locus. Here we consider 
the interpretation of these arc and arclet features. We give a statistical formulation for the 
appearance of the power spectra, based on the stationary phase approximation to the Fresnel- 
Kirchoff integral. We present a simple analytic result for the power-spectrum expected in the 
case of highly elongated images, and a single-integral analytic formulation appropriate to the 
case of axisymmetric images. Our results are illustrated in both the ensemble-average and 
snapshot regimes. Highly anisotropic scattering appears to be an important ingredient in the 
formation of the observed arclets. 
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1 INTRODUCTION 



At low radio frequencies, diffraction from inhomogeneities in the 
ionised component of the InterStellar Medium (ISM) leads to 
multi-path propagation from source to observer. Together with the 
refraction introduced by large-scale inhomogeneities, this leads to 
a variety of observable phenomena — image broadening, image 
wander and flux variations, for example, see Rickett (1990). If the 
source is very small, as in the case of radio pulsars, then the elec- 
tric field has a high level of coherence amongst the various propa- 
gation paths, and persistent, highly visible interference fringes are 
seen in the low frequency radio spectrum. These fringes evolve 
with time, due to the motions of source, observer and ISM. Mostly 
the observed fringes can be understood in terms of wave propa- 
gation through a random medium, with the electron density inho- 
mogeneities being described by a power-law spectrum, close to the 
Kolmogorov spectrum, and physically identified with a turbulent 
cascade (Armstrong, Rickett and Spangler 1995; Cordes, Weisberg 
and Boriakoff 1985; Lee and Jokipii 1976). However, it has been 
known for many years that some pulsars, at some epochs, also ex- 
hibit organised patterns in their dynamic spectra - drifting bands, 
criss-cross patterns and periodic fringes, for example (Gupta, Rick- 
ett and Lyne 1994; Wolsczan and Cordes 1987; Hewish, Wol- 
szczan and Graham 1985; Roberts and Abies 1982; Ewing et al 
1970). The interpretation of these patterns has never been entirely 
clear, although several authors have suggested an association with 
enhanced refraction in the ISM (Rickett, Lyne and Gupta 1997; 
Hewish 1980; Shishov 1974). 

Recently a significant breakthrough in data analysis has oc- 
curred, with the recognition that these organised patterns are more 
readily apprehended by working with the so-called "secondary 
spectrum", i.e. the power spectrum of the dynamic spectrum (Stine- 



bring et al 2001): in this domain the very complex patterns seen 
in the dynamic spectrum typically appear as an excess of power 
along the locus of either a single parabola, or a collection of in- 
verted parabolic arcs (Hill et al 2003; Stinebring et al 2004). Pre- 
vious investigations have demonstrated why parabolic structures 
are generic features of the secondary spectrum (Stinebring et al 
2001; Harmon and Coles 1983), but the properties of the observed 
secondary spectra are not yet understood in detail. In this paper 
we develop a statistical theory of the power distribution in pulsar 
secondary spectra, based on the stationary phase approximation, 
and present some analytic results for simple forms of the scattered 
image. We compare our results with existing data and argue that 
highly anisotropic scattering is required to form some of the ob- 
served secondary spectra. Our results are consistent with, and com- 
plementary to, those of Cordes et al (2003). 

This paper is organised as follows. We begin by developing 
a model secondary spectrum, based on the stationary phase ap- 
proximation, in §2, leading to analytic models for the ensemble- 
average secondary spectra of linear and axisymmetric images (§3), 
and Monte Carlo simulations of snapshot secondary spectra (§4). 
In §5 we consider the appearance of additional, off-centre image 
components; we compare our models with existing data in §6. 



2 THE STATIONARY PHASE APPROXIMATION 

Consider a radio wave propagating over a distance Dp from a 
source to the observer, and passing en route through a phase- 
changing screen at a distance Ds from the observer. We consider 
only the case of a point-like source, implying a large coherent 
patch; this is a sensible approximation for pulsars, which are known 
to be extremely compact sources. Let x denote the 2D position on 
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the phase screen, then the wave ampUtude at a point r in the ob- 
server's plane is 



'(i(r) = 



2-Krl 



d X exp(i$), 



and to high accuracy the phase, $, is given by 



$ = 0(x) + 



(x - /3r) 
2rl 



(1) 



(2) 



with rl := /3Ds/k, = 1 - Ds/Dp, and k = 2-kv/c for fre- 
quency v. Here the pulsar is assumed to lie at the origin of the trans- 
verse coordinates (r and x). In the trivial case where = every- 
where this yields unit wave amplitude at the observer — the usual 
spherical wave amplitude variation, exp(i kz) / z, having been ab- 
sorbed into this normalisation. 

For many purposes the exact integral formulation of eq. 1 can 
be adquately represented by a sum over a finite number of discrete 
points - the stationary phase points - for which the derivative of 
the exponent with respect to the variable of integration vanishes 
(see, e.g., Gwinn et al 1998). This replacement can be made be- 
cause these points, and their immediate surroundings, dominate the 
integral: where the phase is not stationary, the integrand oscillates, 
changing sign rapidly, and there is little net contribution to the in- 
tegral. The condition for a point of stationary phase is V$ = 0, or 



x-/3r 



0. 



(3) 



For each solution of equation 3, x = Xi(r) with i = 1...N, we ex- 
pand the integrand to second order in x about the stationary point. It 
is assumed that these points are all well separated, so that this pro- 
cedure may be applied to each point separately. The integral taken 
around the i th point gives a contribution 



Mi(r) = exp(i$i), 



(4) 



where /i; is the "magnification", and $i the phase for each path. 
The magnification is determined from the phase curvature intro- 
duced by the screen: 

= (1 - -7^, (5) 
in terms of the convergence, /t, and shear 7 := 7I, with 



and 



71 = TT 



72 



dxdy' 



where x = {x, y). The total electric field is then 



(6) 



(7) 



(8) 



and the total intensity is 



I{i>-,t) = u'u '. 



^/Jlijij cos ^ij, (9) 



with $ij s $i — $j. This result is valid not only for the case 
of strong scattering, where N 1 - with each of the A'^ points 
coinciding with a speckle in the image - but also for lens-like (re- 
fractive) behaviour where there may be only a very small number 
of stationary phase points. 



The result we have just derived is a description of the dynamic 
spectrum of the source, because each of the <I>i is a function of 
frequency and time. We will assume that the magnifications change 
only slowly in comparison with the phases. In some circumstances 
- when the point under consideration is close to a critical curve 
(corresponding to the source lying close to a caustic) - the /Xi are 
expected to change rapidly with time and frequency; we will not 
consider such cases, and henceforth we take the /ij to be constant. 
We can approximate the 3> variations as linear in frequency and 
time, over a total observing time At, centred on to, and bandwidth 
Av centred on vq: 

- + ^(t - to) + ^{t^ - uo), (10) 

where = to). Taking the Fourier Transform of equa- 

tion 9 then yields / := I{fu, ft)'- 



At Av 



V/u»ft{exp[i*ij] sinc[7rAt(/t + Ft,ij)] 



X s\nc[K Au{f^ + F^^ij)] + exp[-i<l>jj] 

X sinc[7r At (/« - sinc[7r A;/ (/, - F,,i,)]}. (11) 

For large At, Av this result can be approximated by 

1 ^ 



+ cM-i-^l] S{ft - Ft,,,) 5(U - (12) 
where 6 denotes the Dirac Delta Function, and 



1 



■^1^,1 j — 



1 S*, 



(13) 



27r dt ' 27r df ' 

Our model "secondary spectrum" is therefore (the power spectrum 
of equation 9): 



P{f., ft) = 1*1 ^ E A'iW [S {ft + Ft,ii) 



(47r) 



X <5 (/. F^,ij )+5{ft-Ft,ii) 5{f^-F^,ij)]. 



(14) 



The foregoing analysis shows that the various interference 
terms contribute power to the secondary spectrum in proportion 
to jJLijJLj, at specific fringe frequencies given by the derivatives of 
the phase differences, We note that the secondary spectrum is 
symmetric under the operation {fv, ft) {^fvy~ft) - This arises 
because livyt) is a real quantity, reflecting the interchange asym- 
metry = —^ji. Another manifestation of the synnmetry is the 
fact that fringes with wavevector {fu, ft), and those with wavevec- 
tor {—fv, — ft) are indistinguishable in the dynamic spectrum. Be- 
cause of the symmetry in P, we will henceforth consider only the 
half-plane fi, ^ 0, for which the secondary spectrum just derived 
can be rewritten as 

JV 

P{f^, ft) oc E miij S {ft - Ft,ij) 5 (/. - F.,i,) , (15) 
and this is the form we shall make use of subsequently. 



2.1 Phase relationships in multipath propagation 

The trajectory associated with the z th stationary phase point ex- 
hibits a phase $j = $(xj, r, z/, t), given by equation 2, relative to 
the direct line-of-sight to the source in the absence of a screen. The 
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derivatives of <& with respect to these four variables can be eval- 
uated as follows. Quite generally, differentiating equation 2 with 
respect to r yields 

'f) '''' 

Similarly, differentiating equation 2 with respect to z/ gives 



(x - Prf 



2rl 



2<j) 



(17) 



where we have made use of the result 

'd log 



d log vj^ 

which is appropriate for radio-wave propagation in the Galaxy 
where the refractive index is dominated by the free-electron con- 
tribution. Differentiating equation 2 with respect to x gives 



= 0, 



(19) 



where the final equality (equation 3) applies for stationary phase 
points. Finally, differentiating equation 2 with respect to t yields 

dt L...^ \dt 



(20) 



We will assume that the phase structure in the screen is "frozen", 
and thus is convected with the screen velocity, v^, so that 







(21) 



for all paths. Using the stationary phase condition then gives us 



(22) 



Making use of equations 16 and 22, it is a straightforward ex- 
ercise to evaluate the temporal fringe frequency for a given station- 
ary phase path: 



dt 



dt 



dr 



= — (x-/3r)- (vl-/3vo), 



(23) 



where the observer's velocity is vj,. If the source velocity is Vp, rel- 
ative to some reference frame, then the actual space velocities of the 
screen and observer are ■Vs,o = v',.o + Vp. Only the components of 
these velocities perpendicular to the line-of-sight are of relevance 
here. Defining vx s — /3vo — (1 — /3)vp, where all velocities 
are to be understood as being transverse components, allows us to 
write (dropping the subscript i which denotes the stationary phase 
path under consideration) 



(24) 



where A = 27r//c, and© = (x — I3r)/Ds is the angular separation 
between the path under consideration and the direct line-of-sight to 

the source. 

Implicit in this result for the temporal fringe frequency are the 
Doppler shifts due to the motions of source, screen and observer, 
each of which contributes to v^. All of these terms are linear in 
the relevant velocity and in the angular separation of the paths un- 
der consideration, so for speeds ~ 300kms~^, typical of pulsars, 
and angles of order a couple of milU-arcseconds, the fractional fre- 
quency shifts are ~ 6v/c ~ lO"^'^, hence beat periods of order 



two minutes at a radio frequency of 1 GHz. We note that these 
frequency shifts have not been calculated from a relativistic formu- 
lation, and are therefore not correct to second order in v/c. 
The spectral fringe rate for stationary phase paths. 



27r \'dv 



(25) 



follows directly from equation 17. Again dropping the subscript 
which denotes the particular stationary phase path under consider- 
ation, we have 



2cj3 



(26) 



The physical interpretation of the various terms contributing to F^, 
is straightforward. The first term in equation 26 is the geometric 
contribution to the delay along the path under consideration, while 
the second term is the delay due to the propagation speed of the 
radio-wave within the phase screen; the former is frequency inde- 
pendent, while the latter is dispersive. The geometric delay is ex- 
pected to be ~ d^D/c 10^'' s, for refraction angles 6^2 milli- 
arcseconds, and source distance ~ 1 kpc. 

Equations 24 and 26 apply to the individual stationary phase 
paths, from which wecan compute Ft, ij, F„,ij, via = ^i — ^j. 
Having identified the physical interpretation of the terms contribut- 
ing to Ft, Fv, as being Doppler shifts and delays, respectively, 
we see that the secondary spectrum can be regarded as a "delay- 
Doppler diagram", with the power appearing at locations corre- 
sponding to differences in delay/Doppler-shift between the various 
interfering paths (Harmon and Coles 1983; Cordes et al 2004). 



2.2 Neglect of dispersive delays 

From equations 24 and 26 it can be seen that if the dispersive delay 
(the last term in eq. 26) is neglected, then the fringe frequencies 
obey a quadratic relationship because Ft depends linearly on d, 
whereas F^ varies quadratically. Neglect of the dispersive delay is 
a good approximation, for stationary phase paths, if 



(x-/3r). 



V(/> 



»4. 



(27) 



In other words the screen phase, 0, must be changing on a length- 
scale (L) which is small compared with the distance of the path 
from the optical axis (the direct line-of-sight to the source). The 
meaning of this result can be seen most easily if we consider two 
limiting cases: purely scattering (diffractive), and purely lens-like 
(refractive) phase screens. In the pure scattering case, V0 changes 
sign on a length scale ^ L <^ rp much smaller than the typical 
separation |x — /3r| > rp- Consequently the neglect of dispersive 
delays is a natural approximation to make in the case of a pure 
scattering screen. However, the converse is not necessarily true; to 
see this we can consider a screen in which the phase increases as 
(x — /3r)". In this case the condition given in equation 27 becomes 
n ^ 4, which clearly can be satisfied even though there is no small- 
scale phase structure (no scattering) in the phase screen. 

In the major part of the present paper we will confine our- 
selves to consideration of the case where the dispersive delays are 
negligible, returning to the issue only in §5. This restriction is mo- 
tivated in large part by the attendant simplification of results, and 
encouraged by the fact that the resulting model seems to offer a 
fair description of many of the observed phenomena. Noting that 
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the fringe frequencies Ft.ij and Fu^j include common normalis- 
ing factors (eqs 24 and 26), it is useful to introduce the following 
coordinates (Stinebring et al 2001; Hill et al 2003): 

<l-=r^ft = 01.-62., (28) 
|vx| 

and 

p..= ^f^=el-el (29) 

For simplicity we have chosen the a;-axis to lie along vx. 



3 STATISTICAL MODEL OF THE SECONDARY 
SPECTRUM 

Distant sources (D >1 kpc) observed at low radio frequencies are 
expected to be in the regime of strong scattering, for which there 
are many paths of stationary phase from source to observer. In this 
circumstance it is appropriate to employ a statistical treatment of 
the power distribution in the secondary spectrum. In this section 
we present such a treatment, with emphasis on the analytic results 
which can be obtained for some specific cases. These results should 
prove to be useful for observers undertaking quantitative analysis 
of the power distribution in secondary spectra. 

The intensity distribution in the scattered image can be writ- 
ten in terms of a probability distribution, g(0), for the power in 
the stationary phase paths. Equation 15 shows that the various in- 
terference terms contribute power to the secondary spectrum in 
proportion to /ii/ij, and the expectation value of this quantity is 
g{6i)g{62) (f 61(^62, where we have replaced the specific labels 
i,j hy the generic pair label 1,2. Including only interference terms 
which contribute to the particular fringe frequencies p, q consid- 
ered, the expected power distribution in the secondary spectrum 
can therefore be described by 

P = J <f 61(^62 <?(6'i)g(6>2) S{p - el +el)5{q - 01. + 02.). 

(30) 

Strictly this formulation should be interpreted as yielding the 
ensemble-average (i.e. long-term time-average) of P. However, in 
order to reach this regime observationally one must average data 
over time-scales which are long in comparison to the refractive 
time-scale, and this is generally not the case in practice. In prin- 
ciple one could average the secondary spectra for a given target 
observed over a period of several months or years, much longer 
than the refractive time-scale, but in practice this is not a sensible 
procedure because the large-scale distribution of power in the sec- 
ondary spectrum can undergo profotmd changes on time-scales of 
weeks, indicating that the statistical properties are not stationary 
on these time-scales, so that the ensemble-average regime cannot 
be reached. Instead we are obliged to deal with data in the limit of 
almost instantaneous sampling, or else in the average-image regime 
with only a short averaging time. The various regimes are charac- 
terised by an averaging time, t, which stands in the following rela- 
tionships to the diffractive time-scale, td, and the refractive time- 
scale, tr (Goodman and Narayan 1989): t ^ tr, ensemble aver- 
age; tr ^ t > td, average; t < ta, snapshot. For pulsars observed 
at frequencies v ~ 400 — 1000 MHz, representative time-scales 
are td ^ minutes, tr ^ weeks; a single observation might last for 
t ~ 30 minutes, and is therefore generally in the average-image 
regime. Some data may be in the snapshot regime, and this regime 
is considered in §4. The expected image-plane power distributions 
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Figure 1. "Secondary spectrum" (or "delay-Doppler" diagram) for a linear 
image, at an angle = to the a;-axis (defined by the effective transverse 
velocity vector), with a Gaussian intensity profile (defined by equation 39). 
The axes are expressed in units of Oo = 1/kSo, for q, and 9^ for p. Contour 
levels are set at intervals of 3 dB, and the lowest (outermost) contour is at 
-60 dB relative to P(0.1, 0.1). 

in the average and ensemble-average regimes are identical, and thus 
the formulation given in equation 30 is appropriate for much of the 
data. However, when comparing our models with data we must bear 
in mind that some deviations are expected simply because the data 
are not ensemble averages. 

3.1 Image-plane probability distributions 

In order to compute the secondary spectrum via equation 30, we 
need to specify the probability distribution g(^). The expected an- 
gular distribution of power is given by the spatial Fourier Transform 
of the mutual coherence function (Lee and lokipii 1975): 

g{0) oc j <fp exp{-ikp ■ 6) Tiip). (31) 

For scattering by turbulence in a single screen, the mutual coher- 
ence function is 

V2{p)=eM-D{p)), (32) 

where the phase structure function is given by 

D(p) = ([,^(x + p)-0(x)f). (33) 

In the present paper we will consider only structure functions which 
are power-law in the separation p, and we parameterise these forms 
in terms of the field coherence scale, So'- D{so) = 1. A character- 
istic angular scale for the image is then do = l/(fe So). 

3.2 Highly anisotropic images 

A useful model for highly anisotropic images follows when one 
considers the limiting case of a linear image. A linear source mak- 
ing an angle ip with respect to the x-axis corresponds to an image 
with 

g{6) = h{9x cos tp + Oy sin ip) 6(—9x sin ip + 6y cos ip), (34) 

where the function h denotes the variation in intensity along the 
line. On inserting eq. 34 into eq. 30, there are four integrals and 
four 5-functions. The probabiUty P{p, q) is then found by solving 
the four simultaneous equations, impUed by the four J-functions, 
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Figure 2. Secondary spectrum for a linear image, as figure 1 but for tlie 
image intensity profile given by equation 4 1 (representing the case of scat- 
tering by highly anisotropic Kolmogorov turbulence). Note the difference 
in the scale of the axes relative to figure 1 ; the contour levels are identical 
to those in figure 1. 
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Figure 3. "Secondary spectrum" (or "delay-Doppler" diagram) for an ax- 
isymmetric image with a Gaussian intensity profile (defined by equation 
47). The axes are expressed in units of 9o = 1/kso, for q, and 9^ for 
p. Contour levels are set at intervals of 3 dB, and the lowest (outermost) 
contour is at -60 dB relative to P(0.1, 0.1). 



for 6± — 9x± cos tp + 9y± sin ii and identifying P(p, q) as 
h{9+)h{6-) divided by the Jacobian of the transformation to the 
variables specified by the J-functions. The four simultaneous equa- 
tions give 

_p cos^ ^ + 



and hence 



2q 

2 / 2 

p COS tp — q 
2'q ' 



e± 



6y+ = tanil>Ox+, 
6y- = tanipOx-, 

p cos^ V' =t 9^ 



2q cos ip 

The Jacobian is hence one finds 

h{e+) h{e- 



P{j>,q) = 
with 6± given by eq. 36. 



kl 



(35) 



(36) 



(37) 



3.2.1 Specific examples of anisotropic images 

Here we take the general formulation given above and apply it to 
two specific cases. For simplicity we calculate only the case tp = 
0, corresponding to an image which is elongated along the same 
direction as the effective velocity v_l. 

In the case of a quadratic phase structure function. 



(38) 



with no dependence on py, the resulting image is Gaussian, with 
1 



h{e) = 



27r6>< 



■exp(-eV4e^). 



(39) 



The secondary spectrum corresponding to this image profile is 
shown in figure 1, with q measured in units of 9o, and p in units 
of Oi. The symmetry of the secondary spectrum is due to the sym- 
metry of the image around Ox ~ 0. The contours in figure 1 are 
spaced at intervals of 3 dB, with the lowest (outermost) contour 
corresponding to — 60 dB relative to the peak — the features shown 
here are very faint. In this, and subsequent figures, we have taken 
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Figure 4. Secondary spectrum for an axisymmetric image, as figure 3 but 
for the image intensity profile given by equation 48 (representing the case 
of scattering by isotropic Kolmogorov turbulence). Note the difference in 
the scale of the axes relative to figure 3; the contour levels are identical to 
those in figure 3. 



the peak value of the secondary spectrum to be P( 0.1, 0.1); nor- 
maUsing to these coordinates is somewhat arbitrary, but normalis- 
ing at (0, 0) is not sensible because P — » oo there. The fact that P 
peaks for small p, q simply reflects the fact that the sum of the self- 
interference terms is larger than the sum of cross-power terms in 
equation 15. The most striking feature of figure 1 is the absence of 
power in the region g = 0, p > 0; this is a consequence of the fact 
that we are considering a linear image. The case q = corresponds 
to interference between points which have the same value of 9x, but 
we are considering an image where all of the power is concentrated 
along a line with 9y (x dx, so the secondary spectrum is devoid of 
cross-power at small q — there is only the self-interference which 
appears at p = 0. 

A particularly interesting case to consider is that where the 
structure in the phase screen arises from highly anisotropic Kol- 
mogorov turbulence, in which case 



5/3 



(40) 
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(Narayan and Hubbard 1988), and the image profile cannot be ex- 
pressed in a simple analytic form. Numerically we find that in this 
case the intensity profile assumes a power-law form at large angles, 
with a power-law index of —8/3, and we adopt the approximate 
form 



h{0) 



h{Q) 



1 + (e/eo)8/3 ' 



(41) 



where /i(0) = 4 sin(37r/8)/37r. This image profile yields much 
more power scattered to large angles than the Gaussian profile con- 
sidered above. The corresponding secondary spectrum is displayed 
in figure 2; note the large increase in the area plotted, relative to 
figure 1, even though the contour levels are identical in the two 
figures. This feature is a direct consequence of the fact that the 
Kolmogorov phase screen yields an image with much more power 
at large angles than the Gaussian profile, yielding much larger de- 
lays (i.e. much larger p) at a given contour level. Figure 2 displays 
the same absence of power in the region p > 0, g = as seen in 
figure 1, and for the same reasons. A new feature in figure 2, how- 
ever, is the parabolic form p — delineated by the contouring at 
p,q^ 1. This can be understood by recognising that the dominant 
contributions to secondary spectrum are due to interference terms 
which include the highest intensity region of the image, namely 
1^1 <SC 9o in this case. Referring to equation 36, we see that this 
corresponds to \0- \ <C Oo, hence \p — < 2\q\. More generally, 
for a linear image making an angle ijj with the a;-axis, at large val- 
ues of p the secondary spectrum will exhibit contours which cluster 
around p = sec^ ip. 



3.3 Axisymmetric images 

A priori perhaps the most plausible model is one in which the prob- 
ability does not depend on direction: g(6) = g{9), i.e. the axisym- 
metric case. In this case two of the integrals in eq. 30 can be per- 
formed over the (5-functions, and a useful result requires that one 
perform a third integral analytically. This would then leave a sin- 
gle integral that can be evaluated either numerically, or with the 
use of simple analytic approximations. In view of the symmetry, 
the choice of polar coordinates is the most appropriate. With this 
choice eq. 30 becomes 



P(p,q) = J cieiOidiJi J (16262(14,2 g(ei)g(e2) 

X S{p - el + el) S(q - 01 cosV'i + O2 cosV'2). 
The integrals are carried out in the Appendix; the result is 
(16292 g{e2)g{^/JT^ 



(42) 



P{p,q) 



'■I 

Jo 



xK 



4VPT^6»2 




de2 92 g{62)g{^/^+el) 

[^/¥+^02Y/^ 



l/2\ 



4^p + e^ 02 



1/2N 



(43) 



where K is the Complete Elliptic Integral, and the upper limit of 
integration is = (p — g^)/4|g| + \ {p — q^)/4q\. This result is 
not amenable to further analytic manipulation, in general, and is 
our final form for general axisymmetric images. In some cases it is 
possible to approximate equation 43 with a simpler analytic result. 



For example: if the probability distribution, g{0), is very sharply 
peaked at 6 = 0, then for p,q > and p > q^ the second integral 
in eq. 43 can be neglected (see equations A. 11 and A. 12), and we 
can make the approximation K ~ 7r/2, so that 



'■f 

Jo 



d02 02 g{e2) g{y^) n 1 giVP) 



Vp-q'^ 



\/v-q^ 



(44) 



However, in general the model secondary spectrum for any given 
axisynmietric image must be evaluated numerically from the for- 
mulation given in equation 43 



3.3.1 Specific examples of axisymmetric images 

For the general case of an isotropic, power-law structure function: 

(45) 



we have 



g{d) oc j d^p exp[-ikp-e- D{p)], 

roo 

oc / dttJoi0t/0o) expi-t"). (46) 

Jo 

For Q = 2 it is straightforward to show that the resulting image 
profile is Gaussian, with 



9(0) 



1 



4n0': 



■ exp(-eV4eo)- 



(47) 



The secondary spectrum corresponding to this image is shown in 
figure 3; comparing this to figure 1 (the result for a linear, Gaussian 
image) we see that the structure is broadly similar in the region 
p < g , but the axisymmetric image exhibits much more power in 
the region p > q^. This can be understood in the following way: 
whereas the linear images considered in §3.2.1 possess a unique 
relationship of the form 0y oc 0x, so that a small value of q neces- 
sarily corresponds to a small value of p, no such relationship exists 
for the axisymmetric image — even for q = (so 0xi = 0x2), 
there are many contributing values of 0y2, and hence there is 
power spread over a range of values of p. 

By contrast to the case a = 2, for isotropic Kolmogorov tur- 
bulence, with a — 5/3, the integral does not yield a simple ana- 
lytic form. As with the highly anisotropic case in §3.2.1, evaluat- 
ing numerically reveals a power-law variation in the image surface- 
brightness, with g{0) oc 0~^^^^ for ^ 0o — there is much more 
power scattered to large angles in the case of a Kolmogorov spec- 
trum of turbulence than for the case of a quadratic phase structure 
function. We will use the following approximation for the case of 
isotropic Kolmogorov turbulence: 



g{0) 



l + (6'/6io)ii/3' 



(48) 



where g(0) = 11 sin(67r/ll)/67r'^. The corresponding secondary 
spectrum is shown in figure 4; note the much expanded scale on 
both axes, relative to figure 3, and the enhanced power around the 
parabolic arc defined by p = q^. Comparing to the case of highly 
anisotropic Kolmogorov turbulence, shown in figure 2, we find that 
the axisymmetric case is similar to the anisotropic case forp < q^, 
but exhibits much more power at p > q^. This same feature was 
fotmd in the comparison between linear and axisymmetric Gaus- 
sian images (see above), and occurs for the same reason. 



Parabolic arcs 1 



30 

25 : 

20 : 

15 : 

10 : 

5 : 
_ 



-5 





q 



Figure 5. Secondary spectrum (delay-Doppler diagram) for a linear Gaus- 
sian image in the speclile regime (with i/i = 0, so the image is elongated 
along the direction of vx )■ As for the earlier figures, the axes here are given 
in units of 60 and 9^ , for q and p, respectively. The expected image profile 
is given by equation 39, and the actual image profile is shown in the inset 
at the top of this figure. The ensemble-average secondary spectrum for this 
circumstance is shown in figure 1. 
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Figure 6. As figure 5, but for the case of an axisymmetric Gaussian im- 
age (as per equation 47) in the speckle regime; the actual speckle image is 
shown in the inset box. The ensemble-average secondary spectrum for this 
circumstance is shown in figure 3. 



4 SNAPSHOT SECONDARY SPECTRA 

The results described in §3 are expectation values for the secondary 
spectrum in various circumstances; they represent the ensemble av- 
erage - i.e. the average over many realisations of individual phase 
screens - and are thus appropriate to the long-term time-averaged 
value of observations of secondary spectra. As discussed in §3, in- 
dividual observations of a real source are not expected to be in this 
regime. Usually the results which are reported lie in the average- 
image regime - for which the model developed in §3 is still appro- 
priate, albeit with the caveats given there - but in some instances 
the data may be in the snapshot regime. In this case there is only a 
single realisation of the randomness in the phase screen contribut- 
ing to the observed secondary spectrum, and the phase screen can 
be represented by a particular set of stationary phase points whose 
specific locations are unknown. Obviously the snapshot secondary 
spectra cannot be predicted in detail, but it is important to under- 
stand their general appearance. We give two examples of snapshot 



secondary spectra, for the cases of linear and axisymmetric Gaus- 
sian images; these results may be compared directly with those of 
§3.2.1 and §3.3.1, respectively. 

4.1 Monte Carlo realisations of snapshot spectra 

To demonstrate the appearance of a secondary spectrum from a 
single realisation of a phase screen, rather than the long-term av- 
erage value, we return to the formulation given in §2. To generate 
a secondary spectrum from equation 15 we need to specify the lo- 
cations and magnifications for each of the stationary phase points 
contributing to the image. If the expected image profile is specified, 
e.g. a Gaussian profile in the case of a quadratic phase structure 
function, then we can generate an individual realisation, consistent 
with the expected properties, by the Monte Carlo method. The ex- 
pected magnification of a stationary phase point is independent of 
its location for the models we are employing in this paper, and for 
simplicity we have adopted unit magnification for all the station- 
ary phase points. The coordinates of the stationary phase points 
can be generated as Normally-distributed random numbers, using 
a commercial software package, in either one or two dimensions. 
The results, for an assumed number of N = 100 contributing sta- 
tionary phase paths, are shown in figures 5 and 6 for one- and two- 
dimensional Gaussian images, respectively. 

Figures 5 and 6 show a much richer structure than is evident in 
their ensemble-average counterparts, figures 1 and 3, respectively, 
although the overall power distribution in the snapshots appears 
broadly consistent with that in the ensemble average cases — as it 
should, of course. Of particular interest is that both snapshots man- 
ifest inverted parabolic features; we shall refer to these features as 
"arclets". The arclets can be understood in the following way. For 
the linear image, the speckles are characterised by their locations 
on the X-axis: 9i, i = 1...N, with the interference appearing in the 
secondary spectrum at qij = 9i — 9j,Pij = 9f — Oj.If we consider 
a fixed value of i and allow j to vary, then we recognise that the 
largest value of p corresponds to j — k such that \9k\ < \0j\ for 
all j 7^ k, so that 9k lies close to the origin. For this point we have 
Qik ~ 9i,pik ~ 9f, lying close to the apex of the inverted parabolic 
arclet. Relative to this point, we can determine the locations of the 
other interference terms, by noting that 



qij — qik = — 6 



Hj — Pik = 



whence 



Pa — Pik = —{qij — qik) 



(49) 



(50) 



which describes the locus of points in the ith arclet as the index j 
is allowed to vary. Alternatively, if we vary i and keep j fixed then 
the locus of points corresponds to 



Pij+S'j = {Hj+Ojf, 



(51) 



which is an "upright" parabola of the same curvature (up to a sign) 
as that of the individual arclets. Hence where there are gaps in the 
image - i.e. no stationary phase points around that value of ^ - there 
are corresponding gaps in the intensity of the arclets, and those 
gaps form parabolae cutting through the secondary spectrum. The 
particular case where j = fc in equation 51 describes a locus very 
close to that of the apexes of the arclets: piu ~ 0f ~ so this 
locus also exhibits the same curvature, but passes through the origin 
of the secondary spectrum (t = k). Indeed the self power for all 
image points (j = i), appears at g = 0, p = 0, so all of the 
arclets pass through the origin. These features are strikingly similar 
to what is seen in some of the data (see §6). 
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For the case of a snapshot of an axisymmetric image, the anal- 
ysis leading to equation 50 yields instead the locus 

Pij - Pik = -{qij - qikf - Olj. (52) 

For speckles lying close to the x-axis we have 9yj <C , and these 
points all lie close to the inverted parabolic arclet which defines 
the upper envelope of power for interference with the ith point: 
Pij— Pik ~ —(gij — ga)"^- Referring to figure 6 wc see that there are 
indeed many points lying below each arclet, but the arclets them- 
selves remain clearly visible in many cases. This is because the 
relative delay of the interfering speckles is insensitive to the value 
of 9yj in the region 6yj <C thus leading to a high density of 
power close to the parabolic upper envelope of the arclet. This is 
the same effect, described in §3.3.1, which leads to the visibility of 
the single parabolic arc in figure 4. Similarly, the apexes of the ar- 
clets do not follow a parabolic locus in the case of an axisymmetric 
image; instead they are described by pik = q% + S^j, so p = (f 
forms a lower envelope on the possible apex locations. 
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Figure 7. Secondary spectrum (delay-Doppler diagram) for a multiple- 
component, linear Gaussian image, as described by equation 53. As in fig- 
ures 1-6, q is shown here in units of 6o, and p in units of d^. 



4.2 Evolution of features in tlie snapsliot regime 

The arclets shown in figures 5 and 6 are loci of the interference 
terms between a given stationary phase point (speckle in the im- 
age) and the line 6y =0. (More generally, for a linear image with 
tp ^ 0, the arclets in figure 5 would correspond to interference with 
points close to the locus 6y — 9x tan ip.) In the "frozen screen" ap- 
proximation, the structure of the phase screen does not evolve, and 
consequently individual stationary phase points simply drift across 
the image as the line-of-sight moves relative to the screen (Melrose 
et al 2004). Correspondingly, the apex of each of the arclets seen 
in figures 5 and 6 is expected to drift across the (q, p) plane. In the 
case of an isotropic image the apex of the tth arclet drifts along the 
parabola p = + with \q\ ~ |vj_|/Ds and 9yi = const. 
In the case of a highly anisotropic image making an angle tp = 
with respect to the x-axis, the motion of the apex of the arclets 
is along the parabola p = q^ with |g| ~ \-v±\/Ds. However, if 
ip T^O then the misalignment between the linear image and the ve- 
locity vector means that the stationary phase points cannot remain 
fixed (even approximately) within the screen. Instead, each con- 
tributing stationary phase point must slide in a direction perpendic- 
ular to the image so that at every instant these points constitute a 
linear image. Thus for highly anisotropic scattering we expect the 
apex of each arclet to move along the parabola p = q^ sec^ ip' with 
l^l ~ cos^ -ijj \v±\/Ds. 

Motion of the stationary phase point within the screen, which 
occurs in the case of highly anisotropic scattering with ^ 0, is 
of interest because it implies a short lifetime for the corresponding 
features in the secondary spectrum (even in the "frozen screen" ap- 
proximation). The speed of motion of a given stationary phase point 
relative to the phase structure in the screen is |vx sin i/'l, in a direc- 
tion perpendicular to the image elongation. If the coherence length 
measured in this direction is So - which, by definition, is much 
greater than So, the coherence length measured along the direction 
of image elongation - then the corresponding feature in the sec- 
ondary spectrum is expected to have a lifetime of ~ So/\'v± sm^\, 
and this may be very short unless tp is very close to 0. We note that 
in the development in §3.2, the limit So — » oo was employed for 
simplicity. This is, of course, an idealisation and in practice any 
So ^ So suffices to create highly anisotropic scattering. 

The short Ufetime of stationary phase points in highly 
anisotropic scattering with ip ^ Q suggests that there may be an 



observational bias favouring scattering screens with ■)/) ~ 0, be- 
cause these are the screens which yield persistent high-definition 
arclets. (There may also be further observational biases favouring 
^ ~ in the case of scattering screens of finite spatial extent, be- 
cause of the longer interval over which such screens can contribute 
to the received radiation.) 

The foregoing considerations relate to the evolution of fea- 
tures in the secondary spectrum under the assumption of a phase 
screen which is "frozen", i.e. whose phase profile has no explicit 
temporal evolution. However, the frozen approximation may be 
poor over the time-scale taken for a given point to move right across 
the image. If the phase profile of the screen evolves as a result of in- 
ternal motions or wave-modes with velocity dispersion a, then the 
structure of a coherent patch, of size So, is expected to evolve on a 
timescale ^ So/cr and this gives us an upper Umit on the lifetime 
of any feature arising from a single stationary phase point. 



5 DISTINCT IMAGE COMPONENTS 

To model cases where there are multiple image components, dis- 
tinct from the main image centred onO = (0, 0), the anisotropic 
image model of §3.2 is suitable. Here we consider the simplest mul- 
tiple component configuration, in which there is only one additional 
image, and this component has the same scattering properties (same 
characteristic width, 9o) as the main image; we have chosen an im- 
age with tp = and 

99 1 

ioo''(^^ + Too''(^"^°^'')' ^^^^ 

where h is the Gaussian profile given in eq. 39. In this case the ad- 
ditional image component is weak (1% of the main component) and 
centred on 9x = W9o; figure 7 shows the secondary spectrum for 
this configuration. Compared with figure 1 , the secondary spectrum 
for a single-component linear Gaussian image, there are two new 
features in figure 7: the thin linear feature emerging from the origin, 
and the elongated blob of power with a peak at q = 10, p = 100. 
The presence of the latter peak is expected, based on the analysis 
presented for speckles in §4, as it corresponds to the central portion 
of the weak component interfering with the central portion of the 
main component. However, in contrast to the thin arclets present in 
the snapshot example of figure 6, we see that the image represented 
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by eq. 53 yields a broad swath of extra power extending both to 
larger and smaller p,q. 

We can understand figure 7 in terms of the snapshot proper- 
ties (figure 6) by considering the additional image component to 
be made up of a large number of speckles, each of which yields a 
thin arclet but with the apex at a different location for each speckle; 
it then becomes clear that it is the angular extent of the additional 
image which causes the power to be smeared out in the secondary 
spectrum. Recall that the angular extent was chosen to be the same 
as that of the main image component (eqs. 39, 53), so the main im- 
age and the additional component should exhibit the same spread 
in Doppler shift (q), but the quadratic mapping between image co- 
ordinates and delay (p) means that the additional image component 
exhibits a very extended delay spread because it is not centred on 
the origin. It is helpful to illustrate this numerically, considering 
only the central region of each image component: for the main im- 
age component, extending over Ox = {Oi: l)0o, the delay (relative 
to (0, 0)) ranges over to 1, whereas for the additional image com- 
ponent 9x = (10 ± l)9o with a delay range of 80 to 120. This 
large spread is directly responsible for the large extent in p of the 
additional power in figure 7. 

The fact that the delay varies linearly across the off-centre im- 
age component means that its self-interference term differs from 
that of the main component, even though the two components have 
the same shape. At the centre of the additional image component, 
we have dp/dq = 2q = 20, so that the self-interference of this im- 
age component resides in a linear feature sprouting from the origin 
of the secondary spectrum, close to the line p = 20q . 

Although real data do sometimes show an extended swath of 
power in the secondary spectrum, spread over a large range of de- 
lays, it is more common to see a number of thin arclets. If these 
arclets are transient (see §4.2) then they may be explicable in terms 
of image speckle (§4), but if not then we need to accommodate 
them within the context of distinct multiple image components. We 
have seen (§4) that highly anisotropic images yield thin arclets, be- 
cause there is a unique relationship between p and q when we con- 
sider interference with a given speckle. However, when we consider 
the addition of an extended, off-centre image component, even in 
the limiting case of a linear image, the arclet will exhibit a thick- 
ness (diameter) Ap ~ 4qA0, where AO is the angular diameter 
of the additional image component. The thickness of the arclet will 
be greater than this if the anisotropy of the core of the image is 
less extreme. The existence of persistent, thin arclets with apexes 
at 9 » 1 would therefore place strong constraints on the model 
we have utilised, requiring that the additional image components 
introduced to explain the arclets are subject to much less angular 
broadening than the core of the image. 



6 DISCUSSION 

How well do our statistical models (§3) match the observations? 
For the dozen bright pulsars studied by Stinebring et al (2004), 
scintillation arcs or related phenomena are detected in all of them. 
This effectively rules out an axisymmetric Gaussian intensity pro- 
file (Figure 3) since it does not have sufficient power in the wings 
of the profile to produce scintillation arcs. (The power distribution 
near the origin of the secondary spectrum can be produced by ei- 
ther a Gaussian or a Kolmogorov profile since they are similar for 
< 00.) A Kolmogorov profile with some degree of anisotropy in 
the image is a much better fit to most of the secondary spectra. For 
example, in Figure 8 we show a secondary spectrum for a pulsar 
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Figure 8. The secondary spectrum of PSR B 1929+10 observed at 321 MHz 
at Arecibo at epoch 2003.65. Contour levels are spaced at 3 dB intervals, as 
for the theoretical results shown in figures 1-4, 7; the lowest level contour 
is ^5 dB below the maximum. The observed frequencies {ft, fv) have 
been converted to (q,p) (as per eqns. 28, 29), and scaled to correspond to 
the axes in figs. 1-7, as described in the text. To do so we made use of 
the measured values Dp = 0.33kpc Vp = 163km for this pulsar, 
together with the assumptions that Vp dominates the contribution to v^ 
and /3 = 0.37, and do — 0.50 mas, which are estimated from the main arc 
curvature (Stinebring et al 2001) and measured decorrelation bandwidth, 
respectively. 



(PSR B1929+10) that consistently shows well-defined scintillation 
arcs with a smooth, symmetric fall-off of power from the origin. 
This appears to match best with some form of Kolmogorov image 
(Figure 2 or 4). Determining the degree of asymmetry in the im- 
age will require a more extensive comparison of model and data 
than we attempt here. A visual comparison of the contour plot with 
Figure 4, however, suggests that this image may be close to being 
axisymmetric. 

The most striking secondary spectra are those with pro- 
noimced substructure, particularly those exhibiting inverted arclets. 
As we discussed in §4, these can be reproduced by a model in which 
point-like image maxima interfere with broader features near the 
origin of the image. In Figure 5, for example, the separation be- 
tween isolated image maxima causes gaps, particularly along the 
delay (p) axis, between the inverted arclets. As shown in Figure 9, 
there is remarkable qualitative agreement between observations and 
the simple point-interference model of §4. The arclets in Figure 9 
are clearly delineated (thin), have apexes that lie along a p oc arc 
that has the same (absolute value of) curvature as the arclets, and 
have one side that passes through the origin. All of these features 
are consistent with our model. Notice that the power asymmetry of 
the arclet (extending further to the right than to the left) is the mir- 
ror image of the power distribution along the main arc, where the 
arc has greater power for —q values, but extends out further on the 
+q side of the plot. This is to be expected if, as in our model, the ar- 
clet is the interference between a point-like feature and the central 
portion of the image on the sky. The observation shown here has 
numerous counterparts for this pulsar over several years of recent 
observations, and there is at least one other pulsar (PSR B 1 133+16) 
that exhibits pronounced arclets at numerous epochs. 

The models in §4 were constructed to describe the effects of 
image speckle in the snapshot regime, and the evolution of the cor- 
responding features in the secondary spectnmi was discussed in 
§4.2. The timescale on which the arclets in B0834+06 change is 
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Figure 9. The secondary spectrum or PSR B0834+06 observed at 321 MHz 
on 2004 Jan 8 with the Arecibo telescope. The greyscale is logarithmic in 
relative power, with white being set 3 dB above the noise floor and black be- 
ing set at 5 dB below maximum power. The arclet pattern seen prominently 
on the right hand side of the plot persisted for more than 25 days during 
2004 Jan and moved systematically up and to the right along the guiding 
parabolic arc shown by the dashed Une. The axis scaling is described in the 
text and can be compared directly with the plots of model data. The axes 
in this plot are scaled in the same way as in figure 8, but using the values 
Dp = 0.72kpc, Vp = 175kms-i, /B = 0.33, and Bo ^ 0.72 mas. 

observed to be very long, in some cases at least. In recent observa- 
tions (Hill et al 2004), a persistent arclet pattern has been seen in 
B0834-I-06 for more than 25 days. (This is much longer than either 
the diffractive timescale, of 2 minutes, or the refractive timescale, 
of 2.6 days, for this pulsar at the observing frequency of 321 MHz.) 
Furthermore, the pattern shifts systematically in the direction from 
— g to +q in a linear fashion (|g| = constant) at a rate consistent 
with the proper motion of the pulsar (vj^ ~ (j3 — l)vp). These 
properties point to a phase profile that is localised in one region of 
the scattering screen and is scanned by the pulsar as it moves across 
the sky; this circumstance is consistent with the evolution expected 
in the case of a "frozen screen", as discussed in §4.2. However, the 
large delays associated with the observed arclets (up to 300 /is) re- 
quire the phase coherence length to be as small as So ~ 10^ cm, 
if they are caused by diffraction. If we assume that the character- 
istic velocity dispersion in the phase screen is u > Cs, the sound 
speed, then with ~ 10 kms~^ we expect that individual sta- 
tionary phase points should evolve on a time-scale < 10^ s. This 
is 1/2000 of (the lower limit on) their observed lifetime, and we 
conclude that the observed features cannot be attributed to isolated, 
individual image speckles, ^despite the superficial similarity of the 
data to figure 5. Instead these features must be attributed to either (i) 
a collection of many stationary phase points, if the corresponding 
image feature has a diffractive origin, or (ii) one or more stationary 
phase points associated with strong refraction. In the former case 
one can imagine that the physical cause might be a localised re- 
gion of strong turbulence, in which the individual stationary phase 
points come and go; the Ufetime of the secondary spectrum fea- 
ture is then determined by the decay of the turbulence and may be 
much longer than that of any individual stationary phase point. On 
the other hand, if the feature is associated with strong refraction 
- i.e. a lens - then the relevant length scale in the phase screen 

^ Strictly we should be considering isolated pairs of stationary phase points 
(Melrose et al 2004), but the essence of the argimient remains the same. 



is ^ rF ^ So and the corresponding timescale for evolution of 
the phase screen is naturally much longer Indeed, in the case of 
a lens the phase structure need not be stochastic and might be in 
a quasi-steady state with no detectable evolution (e.g. Walker and 
Wardlc 1998). In cither case the image component must be tightly 
concentrated in order to avoid smearing the power over a range of 
delays (cf. figure 7). These compact refracting/scattering structures 
may be the same structures which are responsible for the Extreme 
Scattering Events (Fiedler et al 1987, 1994). 

Despite the lack of agreement with the snapshot model, the re- 
sults of §4 are useful in understanding the inverted arclets. A com- 
parison of Figure 9 with Figures 5 and 6 indicates that the image 
must be highly elongated in order to produce the observed arclet 
pattern. This is a general result. No inverted arclet observations 
exist that are consistent with an axisymmetric image. Instead, the 
presence of inverted arclets is always accompanied by an absence 
of power in the q = 0, p > region, and the vertices of the arclets 
lie along or very near a definite p = r/g^ locus, with r/ constant 
over many years. All of these aspects follow naturally from an im- 
age that is highly elongated in one dimension and has numerous 
point-like features that self-interfere to give the observed arclets. 
Furthermore, the constancy of t] indicates that the region of the 
interstellar medium being probed has a consistent directionality as- 
sociated with it, as would be produced by anisotropic scattering in 
the presence of a magnetic field (Higdon 1984, 1986; Goldreich & 
Sridhar 1995). This topic is explored further in Hill et al (2004). 

One other insight emerges from a comparison between the 
models presented here and observations of scintillation arcs (e.g. 
Stinebring et al 2004). Interference between two identically scat- 
tered image components (e.g. Figure 7) is rare. Our modeling 
shows that the interference between two identically-scattered com- 
ponents would spread into a broad swath along the scintillation 
arc due to the wide range of differential delay values present. This 
is rarely seen in existing data, with thin arclets or pointlike inter- 
ference features being much more common. This important item 
needs further exploration. 

It is evident from the development in §§2-4 that secondary 
spectra encode information on the structure of radio images (i.e. 
radio-wave propagation paths) of pulsars. A great variety of struc- 
ture is seen in observed secondary spectra, the vast majority of 
which remains to be decoded. Figures 5, 6 and 9 demonstrate that 
secondary spectra with thin arclets or other pointlike interference 
are particularly rich in information. Referring to equation 15, we 
see that A'' points in the image should yield N{N—l)/2 features in 
the secondary spectrum. Consequently we have 0{N'^) constraints 
but only 0{N) unknowns, so if ^ 1 there is enough infor- 
mation to deduce the image structure from the observed secondary 
spectrum. Each stationary phase point (or point-like structure in 
the scattering screen) yields a measurement of the phase (up to an 
overall additive constant), phase gradient, and the combination of 
second derivatives which determines the magnification (eqs. 5-7); 
all of these quantities being determined at the location x, of the im- 
age point. To the extent that the N image points sample the image 
plane, these measurements therefore allow us to deduce the struc- 
ture of the phase screen, (f>. We emphasise that this analysis does 
not require the assumption that the delays are purely geometric; the 
large number of available constraints should in fact provide power- 
ful tests of any approximations which are used in the inversion. The 
main obstacle is observational: we require secondary spectra that 
show strong arclets or other sharply delineated interference features 
with good signal-to-noise ratio throughout. The data shown in fig- 
ure 9 appear suitable, but we have not yet attempted the inversion. 
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7 CONCLUSIONS 

Approximating tlie Fresnel-Kirciioff Integral by a sum over the sta- 
tionary pliase points of an arbitrary phase screen has allowed us to 
develop a probabaUstic description of the power distribution in pul- 
sar secondary spectra. We have presented a single-integral formula- 
tion of the power distribution for general axisymmetric images, and 
a closed form expression for highly anisotropic (i.e. linear) images. 
These results are descriptions of the (ensemble-)average properties, 
but our development makes it obvious how to generate examples of 
snapshot secondary spectra using the Monte Carlo method, and we 
have illustrated this technique with two examples. Our results clar- 
ify the origins of the parabolic arcs in observed secondary spectra, 
and provide an explanation for the inverted parabolic arclets which 
are also sometimes seen. The observed arclets are due to struc- 
tured, highly elongated images. Long-lived arclets require corre- 
spondingly durable features in the phase-screen, and the data show 
that compact, isolated clumps of refracting/scattering material must 
exist in the interstellar medium. Consideration of the lifetime of 
highly anisotropic image components highlights an observational 
bias in favour of those which are elongated along the direction of 
the effective velocity vector. 
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APPENDIX A: EVALUATION OF INTEGRALS FOR 
AXISYMMETRIC IMAGES 

In evaluating the integral 

p = J <fei<fd2p{di)p{e2)5{p-e\+el)5{q- 91^ + 02.), 

(Al) 

one may assume p without loss of generality, with the func- 
tion for p < determined from that for p > by P{—p, q) = 
P(p, —5). For the axisymetric images under consideration here, the 
secondary spectrum is also symmetric in q, and in this Appendix we 
further limit the domain of our calculation to g 0. One may carry 
out any two of the four integrals in eq. A.l over the (5-functions. 
We evaluate the integrals in polar coordinates, dPOi = OidOidtpi. 
Carrying out the integral over d^9i gives 



(A2) 



where 6iy is determined by the (5-functions, and is given by 

=p+e'2- {q + e2COSij2f. (A3) 

The boundary of the physical region is defined by F = 0. This 

corresponds to a parabola with apex at 92x — {p — q^)/2q, and 
asymptotes at 62y ±00, 62X — > +c». The integral is over the 
region outside this parabola. For {p — q^) /2q > the apex of 
the parabola is to the right of the origin, and the integral is over 
< V2 < Stt for 62 < (j> — q^)/2q, and over cos'!/'2 < cos ?/>+, 
where 



COS'!/'+ 



VP + el 



(A4) 



forS2 > {p—q^)/2q.For {p — q'^)/2q < the apex of the parabola 
is to the left of the origin, so that a region around the origin is not in 
the range of integration, and the integral is over cos ^^2 ^ cos 1/^+ 
for^a > \ip-q^)/2q\. 

One procedure for evaluating the angular integral is to first 
change variables from tp2 to t, using 



t = tan 



1p2 



2dt 



COSV'2 = 



1-e 



Then one finds 



.1/2 



[(t'+a'){t' + b')] 



1/2 



where 



2 

a = 



^/p + ei + q + e2 

\fpT&i -\-q-02 



0'i~q-t 



Vp + 01 -q + 02 



(A5) 



(A6) 



(A7) 



For p > and Q < 62 < {p — q^)/2q (with q ^ 0), we have 
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> b'^ > 0, while for 02 > \ ip - q^)/2q\ we have &^ < < a^. 
Introducing — —h^, and making use of Gradshteyn and Ryzhik 
(1965: 3.152, 8.111.2, 8.112.1, 8.113.1) we find that the requisite 
integrals are 

0- .IkC^^V ,A8) 



Jo 



[(a2+a;2)(a;2 + 62)]i/2 a 

and 

r dx _ 1 

J a [(a2+x2)(x2-d2)]V2 ^^23^2 \,/^^ ) ' 

(A9) 

where K is the Complete Elliptic Integral, with 

K(fc) = f (l + ^ + ...). (AlO) 



Jo 



The foregoing results imply 
for < ^2 < (p - 9^)/2g (with q > 0), and 

for 02^ \{p-q^)/2q\, where 

^= (^y^p + ei+Sa^ -q^ (A13) 

and 



B = 46i2Y'p + ^i- (A14) 

Substituting these results into equation A2 then yields our final re- 
sult for the secondary spectrum corresponding to an axisyimnetric 
image (equation 43). 



